Kondo correlation and spin-flip scattering in spin-dependent transport 
through a quantum dot coupled to ferromagnetic leads 
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We investigate the linear and nonlinear dc transport through an interacting quantum dot con- 
nected to two ferromagnetic electrodes around Kondo regime with spin-flip scattering in the dot. 
Using a slave-boson mean field approach for the Anderson Hamiltonian having finite on-site Coulomb 
repulsion, we find that a spin-flip scattering always depresses the Kondo correlation at arbitrary po- 
larization strength in both parallel and antiparallel alignment of the lead magnetization and that 
it effectively reinforces the tunneling related conductance in the antiparallel configuration. For sys- 
tems deep in the Kondo regime, the zero-bias single Kondo peak in the differential conductance is 
split into two peaks by the intradot spin-flip scattering; while for systems somewhat further from 
the Kondo center, the spin-flip process in the dot may turn the zero-bias anomaly into a three-peak 
structure. 
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Spin-dependent electron transport through systems 
consisting of two ferromagnetic layers (FM ) sandwiched 
by a quantum-dot (QD) has recently become one of the 
major focuses of the rapidly developing spin-electronics. 1 
Different from other nonmagnetic sandwiches, a quantum 
dot features its small size and the strong many-body cor- 
relation among electrons in it such that, not only the 
Coulomb blockade can dominate its electronic transport, 
but particularly, a significant Kondo effect may arise. 
Very rich behavior has been disclosed when varying the 
strength and/or the relative orientation of the spin po- 
larizations of two magnetic electrodes, the dot energy 
level and the on-site Coulomb correlation in the dot. 2-10 
Effect of spin-flip scatterings in the QD has also been 
explored very recently. By means of the equation of 
motion method 11 and the Ng ansatz 12 , Zhang et. al. 9 
predicted that, in a FM/QD/FM system having infinite 
on-site Coulomb repulsion, the presence of the spin-flip 
process in the dot may split the original single Kondo 
resonance in the density of states into two or three well- 
defined peaks. Utilizing a slave boson mean field (SBMF) 
technique for infinte-J7 systems, 13 ' 14 Lopez and Sachez 10 
found a splitting of the zero-bias peak of the nonlinear 
differential conductance when the spin-flip scattering am- 
plitude is of the order of the Kondo temperature. How- 
ever, the failure of detecting the Kondo peak splitting- 
induced by lead magnetization in the parallel configura- 
tion and the unphysical prediction of complete quench- 
ing of the zero-bias anomaly and vanishing widths of the 
nonzero-bias peaks at large spin-flip scattering, 10 indi- 
cate the limitation of the infinite-?/ SBMF approach to 
FM/QD/FM systems. 

In this Letter we investigate effects of intradot 
spin-flip scatterings on the dc conductance through a 
FM/QD/FM system using the finite-?/ SBMF approach 
of Kotliar and Ruckenstein. 15-17 With the help of four 
slave boson parameters, this SBMF approach not only 
allows to deal with finite on-site Coulomb repulsion 



but also to take account of charge fluctuations to cer- 
tain degree. For a FM/QD/FM system without spin- 
flip scattering 8 it derives a clear splitting of the Kondo 
peak in differential conductance due to lead magneti- 
zation in the parallel configuration in agreement with 
other techniques, 6, and predicts the physically reason- 
able peak widths and heights. 

We consider a quantum dot (QD) of single bare level 
ed having a finite on-site Coulomb repulsion U and a 
spin-flip scattering amplitude R, is coupled to two ferro- 
magnetic leads through tunneling Vk aCr (a = L,R). Both 
the left and the right leads are magnetized along the z 
axis but may be in parallel (P) or antiparallel (AP) align- 
ments (configurations) and the electrons in the leads are 
described by wavevector k and spin index a with energy 
£fe aCT (c = ±lor TJ.)- 

In terms of the finite-?/ slave-boson approach, 15 we in- 
troduce four auxiliary boson operators e, p a and d, which 
are associated respectively with the empty, singly occu- 
pied and doubly occupied electron states of the QD, to 
describe the above physical problem without interparticle 
coupling in an enlarged space with constraints. Within 
the mean-field scheme one can start with the following 
effective Hamiltonian: 
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where c\ (cda) and ct a (ck a a) are the creation (anni- 
hilation) operators of electrons in the dot and in the 

electrodes, the three Lagrange multipliers A^ and A^ 2 ' 
take account of the constraints, and in the hopping term 
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the fermion operator c da and c da are replaced by z^c da 
and CdaZa to consider the many body effect on tunnel- 
ing. z a cmbodys the physical process with which an 
electron with spin a in the dot is annihilated: z a = 

(1-dU- vlVa)~^ (eV +pU)(l ~ e f e - vlvsY^ ■ The 
Heisenberg equations of motion for the four slave-boson 
operators derived from the effective Hamiltonian together 
with the three constraints, form the basic equations. 
Then we use the mean-field approximation in the sta- 
tistical expectations of these equations, in which all the 
boson operators are replaced by their expectation values 
which serve as the instrumental parameters in the mean- 
field scheme. In the wide-band limit the final equations 
are: 





E^ + 2A«e = 0, 


(2) 
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J2\P«\ 2 + \e\ 2 + \d\ 2 = l, 

cr 
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J d^G<M = \P<r\ 2 + \d\ 2 , * = ±1. 


(6) 



Here 

K a =± J du{z a Re{G d M)[f L {uj)T LG + f R (u)T Ra }}, 

(7) 

where f a (u>) = l/(e^^ _M °) + 1) is the Fermi function 
and n a is the chemical potential of the ath lead, which 
is assumed in an equilibrium state at temperature 1/(3, 
and T ar7 (w) = 2irJ2k a \Vk a a\ 2 S(w - e ka<7 ) is the cou- 
pling strength function between the QD and the lead a. 

G d ^,(co) and G d<y(yl {uj) are the elements of the 2x2 re- 
tarded (advanced) and correlation Green's function ma- 
trices Q r d a \uj) and G d {uj) in the spin space of the QD. 
In the SBMF scheme, the strong correlating problem of 
the dot electrons with finite Coulomb repulsion U is re- 
duced to an effective noncorrelating one described by 
the effective Hamiltonian (1) with the slave-boson op- 
erators treated as c-numbers. The retarded (advanced) 
Green's function can be expressed in the form renormal- 
ized due to dot-lead couplings: Gj"'y = - J, , in 

° a v ' uil — tl d ±tT 

which H daa = e da (e da = e d + \<r ) and H daal = R 
(<T 7^ cr') take up the dot-level shifting and splitting 
and the spin-flip process, f CT0 ./ = \ \z a ^{V La + Y Ra )5 aa > 
reflects the dot-level broadening, and I is a unit ma- 
trix. The correlation Green's function can be ob- 
tained by the Keldysh equation: G^- = G^E^G^j, with 

The electric current flowing from the left lead into the 



QD is obtained from the rate of change of the electron 
number operator of the left lead: 11,18 

Il = jJ ^E^I^I'Kg^-glJ/l + g^]. (8) 

In the steady transport state, the current flowing from 
the QD to the right lead must be equal to the current 
from the left lead to the QD, and the formula (8) can be 
directly used for calculating the current flowing through 
the lead-dot-lead system under a bias voltage V between 
the two leads: I = II- 

We assume that the left and right leads are made from 
identical materials and that, in the wide band limit, 
the effective coupling strength functions are constants 
for each direction of magnetization, such that in the P 
configuration T R ^(u>) = = Tj and T R i(u>) — 

Tli{oj) = Ti, and in the AP configuration T r ^(lj) — 
T ll (oj) = T x and T RL (ui) = T L1 (u) = T T . We will take 
r = (r^ + T|)/2 as the energy units and define the spin 
polarization as P = (T^ — T|)/(2r). The Kondo temper- 
ature for finite- U system in the case of P = 0, given by 
T°k = *7\//3cxp(-7r//?)/2^ with (5 = -2UT/e d {U + e d ), 
will be used as a reference dynamical energy scale. In 
the following we will study FM/QD/FM systems having 
a fixed finite Coulomb repulsion [7 = 6 and concentrate 
on the effects of the spin-flip scattering on linear and non- 
linear transport. In calculating the dc current under a 
bias voltage V between two leads, we choose the chemical 
potential [i^ = —y> R = eV/2 for the left and right leads. 

The finite-t/ SBMF method was known to yield quali- 
tatively correct conduction behavior at low temperatures 
within the bare dot-level range —1.2U < e d < 0.2(7. 
Fig. 1 shows the calculated zero-temperature linear con- 
ductance G = (dl /dV)v=o as a function of e d of the 
QD having spin-flip scattering of different amplitude R 
at fixed polarization P = 0.7 in AP and in P configura- 
tions. 

In the case of R = 0, the G-vs-e d curve contains three 
regimes, covering the resonance peak due to the dot level 
around e d = 0, the charging peak around e d = —U, and 
the Kondo peak centered at e d = — U/2. In the AP con- 
figuration, all the renormalized parameters are identical 
for up and down spin indices at arbitrary P and electrons 
with up-spin and down-spin are equally available in the 
whole lead-dot-lead system, favoring the formation of the 
Kondo-correlated state within a relatively wide dot-level 
range centered at e d = —3. The G-vs-e d curves appear 
to be gentle-top hump structure. On the other hand, 
in the AP configuration the available minority-spin (e.g. 
up-spin) states in the right lead decline due to P > 
and the transfer of the majority-spin (up-spin) electrons 
from the left lead to the right lead is suppressed by the 
finite polarization, such that the conductance of the sys- 
tem goes down with increasing polarization from P = 
and vanishes at P = 1. The effect of a finite polarization 
is to reduce the height of G-vs-e d hump, while the shape 
of the curve remains essentially unchanged. 
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FIG. 1. The linear conductance G versus e d at P = 0.7 
(left) and versus P at e d = —2 (right) in AP and P config- 
urations for different values of R. The systems have a finite 
value (7 = 6. 

In the P configuration, finite polarization reduces the 
available minority-spin electrons in both electrodes and 
splits the dot level for up and down spins, thus weakens 
the Kondo correlation (narrowing the Kondo peak cen- 
tered at td = —3) and slightly broadens the tunneling 
peaks around = and ed = —U. However, since both 
up- and down-spin electrons can tunnel through from the 
left lead to the right lead, the heights of the Kondo peak 
and tunneling peaks are essentially not affected by the 
finite polarization and the unitary limit (G = 2e 2 /h) can 
still be reached as long as P < 1. The main feature in the 
P configuration is the central Kondo peak progressively 
narrowed with increasing P. 

A coherent spin-flip scattering (R > 0) in the dot plays 
three important roles. (1) It lifts the degenerate dot level 
ed to ed ± R even in the presence of a strong on-site 
Coulomb repulsion, with states being the linear combi- 
nations of the up- or down-spin ones. This level split- 
ting may lead to somewhat broadening of the resonance 
peak and charging peak. (2) It strongly suppresses the 
Kondo correlation, such that the conductance maximum 
at ed = —3, which is fully Kondo-induced, drops sharply 
when R > 0. (3) In the AP configuration, R > makes 
it possible for the majority-spin electrons of the left lead 
to change their spin orientations after tunneling into the 
dot, in favor of their continuing to tunnel into the right 
lead, and thus enhances the tunneling related conduc- 
tance. As can be seen in Fig. 1, the tunneling peaks and 
the charging peaks in the cases of i? = 0.1, 0.2, 0.3 and 0.4 
are all higher than those of R = 0. In the P configuration, 
however, the tunneling-related conductance is essentially 
not affected by the presence of spin-flip process. The G- 
vs-ed curves in Fig. 1 are the result of competition and 
compensation of the above three effects. We see that, 



in the P configuration, the spin-flip scattering reduces 
the central Kondo peak sharply while slightly shifts the 
resonance peak and charging peak apart, finally turning 
the three-peak G-vs-ed curve into a Coulomb blockade 
type double-peak one. In the AP configuration, a deep 
valley develops progressively at the center ed — — 3 from 
the previous G-vs-ed hump, leading also to a Coulomb 
blockade type double peak one with enhanced heights of 
and enlarged distance between two peaks. 

Effects of changing polarization strength on linear con- 
ductance in the presence of spin-flip process are shown 
in the right part of Fig. 1 for the system of = — 2 in 
the Kondo regime. In the P configuration, G always goes 
down with growing P for fixed R and with growing R for 
fixed P, indicating that both the polarization and the 
spin-flip scattering suppress the Kondo resonance. In the 
AP configuration, although the linear conductance G still 
declines with increasing P in the case of R = or with 
increasing R in the case of P = 0, anomaly shows up at 
finite P and finite R due to the competing effects of po- 
larization and spin-flip scattering on Kondo correlation 
and resonance tunneling. At R = 0.1 the conductance 
descends much slower than at R = 0, and G even rises 
with increasing polarization for R > 0.2. 
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FIG. 2. Zero-temperature dl/dV versus eV/T% for 
e d = -2 at P = 0.5 (left) and for e d = -1 at P = 0.7 
(right) having different values of spin-flip amplitude R in AP 
and P configurations. The systems have a finite value U = 6. 



Zero-temperature differential conductance dl/dV is 
calculated at bias voltage eV up to several T^. Fig. 2 
left shows dl/dV versus eV/1% at different spin- flip scat- 
tering amplitudes R = 0,0.1,0.2,0.3 and 0.4 but same 
polarization P = 0.5 for the system of = —2 in the 
Kondo regime (T% = 0.144). For this system at P = 0.5 
in the absence of spin-flip process, the dot-level split- 
ting induced by the lead polarization alone in the P con- 
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figuration has not yet been able to split the broadened 
Kondo peak. We have essentially the conventional zero- 
bias anomaly of Kondo correlation in both AP and P 
configurations for R = 0. Turning on the intradot spin- 
flip scattering suppresses the Kondo effect and lifts the 
dot-level degeneracy of order of R, giving rise to a signif- 
icant drop of the height of the zero-bias anomaly and, up 
from R = 0.1 (~ 0.7T°), dl/dV begins to split into two 
peaks centered around eV — ±2R in both configurations. 
This result is in agreement with the recent prediction of 
Lopez and Sanchez. 10 Note that, however, our dl/dV for 
this system does not show any sign of complete quench- 
ing of the zero-bias conductance and the split peaks pre- 
dicted here exhibit physically reasonable widths of order 
of renormalized level broadening. 

Situation becomes somewhat different for the system 
of Cd = — 1, a little further from the center of the Kondo 
regime, as shown in the right part of Fig. 2. For this 
system at P = 0.7, in the P configuration the lead- 
polarization induced dot-level splitting is already large 
enough to show up as the zero-bias peak splitting in the 
case of R = 0. Introducing spin-flip scattering further re- 
duces the zero-bias differential conductance and enlarges 
the separation between the two split peaks in such a way 
that the distance increment is roughly 4 times of the R 
increment except the case from R = to R = 0.1. Note 
that for this system the zero-bias conductance, which is 
already much lower than the unitary limit 2e 2 /h due to 
the suppression of lead magnetization, comes only partly 
from Kondo correlation and partly from resonance tun- 
neling. The spin-flip process suppresses the Kondo con- 
tribution while keeps the latter essentially unchanged. 
Therefore the drop of the zero-bias dl/dV is modest in 
comparison with that in the left part of Fig. 2. On the 
other hand, the peak splitting develops progressively due 
to the increase of i?-induced level splitting. The appear- 
ance of three-peak structure in the curves of R = 0.2, 0.3 
and 0.4 is just the result of these two effects. In the 
AP configuration, there is no dot-level splitting and the 
differential conductance remains a zero-bias single-peak 
curve for any strength of the lead polarization in the case 
of R = 0. The spin- flip scattering, while suppressing 
the Kondo contribution, helps to enhance the resonance- 
tunneling related conductance. For = — 1 system the 
latter overcompensates the Kondo reduction, such that 
we see the whole dl/dV curve of R = 0.1 rising over 
that of R = 0. Note that at R = 0.1, the spin- flip in- 
duced level splitting has not yet been able to destroy 
the single-peak structure and we have only a broadened 
peak. The R = 0.2 induced level splitting leads to the 
appearance of shoulders on both sides. Two side peaks 
clearly show up at larger spin-flip scattering amplitudes. 
We see three-peak dl/dV curves for R = 0.3 and 0.4 in 
the AP configuration. 

In summary, we have demonstrated that a spin-flip 
scattering in QD always suppresses the Kondo correla- 



tion at arbitrary polarization strength in both P and AP 
alignment of the lead magnetizations, and the effect is 
stronger for P configuration at larger P. It effectively 
reinforces the tunneling related conductance in the AP 
configuration. For systems deep in the Kondo regime, 
the zero-bias single Kondo peak in the differential con- 
ductance is split into two peaks by the intradot spin-flip 
scattering; while for system somewhat further from the 
symmetric point, the spin-flip process in the dot may 
turn the single-peak conductance curve into a three-peak 
form. The predicted positions and widths of the split 
peaks are in rough agreement with the spin-flip induced 
energy splitting and the tunneling induced level broad- 
ening, renormalized by the many-body effect. The main 
physics concerned in this Letter is in the energy scale 
of the Kondo temperature or the effective coupling 
strength T, which are easily accessible experimentally. 
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